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A systematic replica field theory calculations are analysed using the examples of two particular 
one-dimensional "toy" random models with Gaussian disorder. Due to apparent simplicity of the 
model the replica trick calculations can be followed here step by step from the very beginning 
till the very end. In this way it can be easily demonstrated that formally at certain stage of the 
calculations the implementation of the standard replica program is just impossible. On the other 
hand, following the usual "doublethink" traditions of the replica calculations (i.e. closing eyes on 
the fact that certain suggestions used in the calculations contradict to each other) one can easily 
fulfil the programme till the very end to obtain physically sensible result for the entire free energy 
distribution function. 
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I. INTRODUCTION 



In recent years there is a renewed interest to the mathematical status of the replica method widely used in disordered 
systems during last four decades. For the calculation of thermodynamic quantities averaged over disorder parameters 
(e.g. average free energy) the method assumes, first, calculation of the averages of an integer n-th power of the 
partition function Z(n), and second, analytic continuation of this function in the replica parameter n from integer 
to arbitrary non-integer values (and in particular, taking the limit n — > 0). Usually one is facing difficulties at both 
stages of this program. First of all, in realistic disordered systems the calculations of the replica partition function 
Z(n) can be done only using some kind of approximations, and in this case the status of further analytic continuation 
in the replica parameter n becomes rather indefinite since the terms neglected at integer n could become essential 
at non-integer n (in particular the limit n — > 0) (H, HJ. The typical example of such type of trouble is provided by 
the classical Kardar's solution of (1 + 1) directed polymers in random potential where due to the approximation 
used at the first stage of calculations (when the parameter n is still integer) the resulting free energy distribution 
function appears to be not positively defined [!, Q (see also @). On the other hand, even in rare cases when the 
derivation of the replica partition function Z(n) can be done exactly, further analytic continuation to non- integer 
n appears to be ambiguous. The classical example of this situation is provided by the Derrida's Random Energy 
Model (REM) in which the momenta Z(n) growths as exp(n 2 ) at large n, and in this case there are many different 
distributions yielding the same values of Z(n), but providing different values for the average free energy of the system 
Q. Performing "direct" analytic continuation to non- integer n (just assuming that the parameter n in the obtained 
expression for Z(n) can take arbitrary real values), one finds the so called replica symmetric (RS) solution which 
turns out to be correct at high temperatures, but which is apparently wrong (it provides negative entropy) in the low 
temperature (spin-glass) phase. In the case of REM the situation is sufficiently simple because here one can check 
what is right and what is wrong comparing with the available exact solution (which can be derived without replicas). 
Unfortunately in other systems the status of the results obtained by the replica method is much less clear. 

In the case of the mean-field spin-glasses 0| the replica partition function also growths as exp(n 2 ) at large n, and 
its "direct" analytic continuation to non-integer n, as in REM, provides wrong RS solution in the low temperature 
spin-glass phase. Here the solution which is generally believed to be correct is obtained via the Parisi replica symmetry 
breaking (RSB) scheme (in the case of REM it reduces to the special case which is called one-step RSB), and it is 
derived in terms of a heuristic procedure and not as a proper analytic continuation from integer to non integer values 
of n of the replica partition function Q . Recently the results obtained in terms of the RSB scheme has been confirmed 
by independent mathematically rigorous calculations (see Q and references therein) . Although no one seems to doubt 
now that the RSB heuristic procedure provide correct results, the problem is that until now no one was able to explain, 
why it provides correct results? 

Presumably the most notable progress in the studies of the subtleties of the replica method has been achieved 
recently in the context of the random matrix theory, where the remarkable exact relation between replica partition 
functions and Painleve transcendents has been proved [I3 - [l5j . One can also mention here recent exact replica 
solution for one-dimensional directed polymers in random potential, which (unlike all previous examples) did not 
involve an analytic continuation from integer to non-i nteg er replica parameter n and where the results obtained were 
also expressed in terms of the Painleve transcendents HH ■ 

In this paper I would like to consider two examples of the systematic replica field theory calculations using very 
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simple random systems for which, at first sight, every step of the replica program is under control. It turns out, 
however, that even in these extremely simple cases the derivation of the physical results inevitably requires the usual 
replica method "cheating" (the advantage of the simple system is that here one can easily see how it goes). Of course, 
it is not that present research explains something deep about the replica method. The aim of the paper is (once 
again) to turn the attention to the existing paradoxes, and to promote (once again) the idea that "something has to 
be done" : it looks rather uncomfortable that on one hand we have extremely robust method, which in most of the 
cases works perfectly well, while on the other hand, we do not understand why it works. 



The "toy" models considered in this paper are "extracted" from the one-dimensional directed polymers in a quenched 
random potential (for the detailed physical analysis of the obtained results see These systems describe an elastic 

string directed along the s-axis within an interval [0, L] with displacements defined by the scalar field 4>{x) having the 
elastic energy density proportional to (d x (f) 2 . Randomness enters the problem through a disorder potential V[<f)(x), x] 
competing against the elastic energy (see e.g. [I3])- The problem then is defined by the Hamiltonian 

H[4>{x),V\ = J dx{^[d x <f>{x)] 2 + V[<j>(x),x}}; (1) 

The disorder potential V[4>{x),x\ is Gaussian distributed with a zero mean V(cj>,x) = and a correlator 



V(cj), x)V{<f>', x') = S(x - x')U(<f> - (j>'); (2) 

which is defined by a correlation function U(4>). The above equation implies that the random potential correlations 
are "translation invariant" in the <p direction depending only on the difference (4> — </>'). 

General strategy of the replica calculations for this system is in the following. For the string with the zero boundary 
conditions at x = the partition function of a given sample is 

/+oo „ c f > (L)=y 
dy / T>[<f>(x)] e-^< yI (3) 

where the integration goes over all trajectories <fi(x) staring at the origin, and /3 denotes the inverse temperature. On 
the other hand, the partition function is related to the total free energy F[V] via 

Z[V] = exp(-pF[V}) (4) 

The free energy F[V] is defined for a specific realization of the random potential V and thus represent a random 
variable. Let us take the n-th power of both sides of Eq.(j4]) and perform the averaging over the random potential V: 



Z[n,L] = [exp{-PnF[V]) j (5) 
The quantity in the l.h.s of the above equation 

Z[n,L]= (z[V]) n (6) 

is called the replica partition function, and it is defined originally for an arbitrary integer parameter n. Let us suppose 
that at large L the fluctuating free energy of the system scales with the system size as F cx L u , i.e. it is characterized 
by a single universal exponent lu. Redefining F = fL u , we introduce a random quantity / ~ 1 which can be described 
by a distribution function Vl(J) (depending on the system size L). In this way, instead of eq.© we get the following 
general relation between the replica partition function Z[n, L] and the distribution function of the (rescaled) free 
energy fluctuations Vl(/)' 

/+oo 
dfV L (f)e- pnL " f (7) 
-oo 

The above equation is the bilateral Laplace transform of the function Vl(I), and at least formally it allows to restore 
this function in terms of the replica partition function Z[n,L]. In order to do so we have to compute Z[n,L] for an 
arbitrary integer n and then (if the result would permit!) perform analytical continuation of this function from integer 
to arbitrary complex values of n. Introducing a new complex variable 

s = finL" (8) 
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and denoting 



Z[^,L]=Z L (s) (9) 



instead of eq.([7]) we get 

/+oo 
dfv L (f)e- sf (10) 
-oo 

According to this relation the distribution function VlH) can be reconstructed via the inverse Laplace transform 

r+ioo , 

v L {f) = J 2^^( s ) es/ ; (ii) 

where the integration goes over the contour parallel to the imaginary axis. Finally, provided there are exist a finite 
thermodynamic limit function 



lim Z L {s) = Z*(s) (12) 



we can find the distribution function 



r+ioo r 
J-ioo 



(13) 



which would describe the statistics of the rescaled free energy fluctuations / in the infinite system. The above equation 
defining V*(f) contains no parameters and hence is expected to be universal. Note also that according to the relation 
s — finLt* 3 in the thermodynamic limit L — > oo the relevant values of the replica parameter n ~ — > 0, which 
explains why the two limits L — > oo and n — > do not commute [H . 

Usually the problem formulated in eqs.([T])-(G±]) is studied in the context of the short-range correlated disorder 
potential, i.e., for a rapidly decaying function U(4> — > oo) — » (which for simp licity is often replaced by the 5- 
function) In this case the free energy fluctuations scale as L 1//3 0, [l^, [n], Il8l - l20j |. On the other hand, if we would 
like to study the statistics of small string displacements we could develop the random potential V[(f), x] in powers of 
4> <^ 1 keeping the first two terms only: 

V[<j>,x]siVo(x) + Vi(x)<i> (14) 



where Vq(x) and V\{x) are the Gaussian uncorrelated random parameters with the zero mean, Vq(x) = V\(x) = 0, 
and the correlators 



Vo(x)V {x') = vS(x-x') 



Vx{x)V x {x>) = u5{x-x') (15) 



Vi(x)Vo(x') = 

which implies that unlike eq.([2]), the random potential correlator is not translation invariant: 



V{<j>, x)V{cj>',x') = 5{x - x 1 ) [v + ucjxj)'} (16) 
In this way, instead of eq. (JlJ we arrive to the Hamiltonian 

H[4>{x),V}= ( dx{^[d x <t>(x)] 2 + Vi(x)<f>(x)+V (x)y, (17) 



o 

Now, lifting the requirement <f> <C 1 , we are getting simple Gaussian random force model to be studied in this paper 
by the replica method (Section II). 

The thermodynamic limit, L — > oo, of this system has been studied earlier [2l|. Here I would like to concentrate on 
the thechnical details of the calculations. Due to apperent simplicity of the model the replica trick calculations can 
be followed here step by step from the very begining till the very end. In this way it can be easily demonstrated that 
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formally at sertain stage of the calculations the implementation of the replica program (as it is declared at the very 
begining) is just impossible. On the other hand, following the usual "doublethink" traditions of the replica calculations 
in disordered systems (i.e. closing eyes on the fact that certain suggestions used in the calculations contradict to each 
other) one can easily fulfill the programme till the very end to obtain very nice and physically sensible result for the 
free energy distribution function. Moreover, in this particular case we can be sure that the result obtained in this way 
is indeed correct, as it can also be derived via direct calculations without replicas fl6j . It turns out that regardless of 
the apparent simplicity of the model, its free energy distribution function (which will be derived here for an arbitrary 
finite system size L) is rather non-trivial. 



In Section III we consider slightly modified version of the above "toy" system. Namely, instead of developing the 
random potential itself, eq. (fT4| . one can consider the development of its correlation function U ((/>), eq.([2]). Again, 
keeping the first two terms only, one gets 

U{4>) ~ v - X -u4? (18) 

(which, unlike cq. (|16[) , would preserve the translation invariance of the random potential correlations). It turns 
out that this, seemingly rather innocent modification provokes quite dramatic consequences. As before, we lift the 
requirement <fi -C 1 and accept the above "truncated" correlation function, eq. (fT8|) . as valid for the whole range of the 
scalar fields <f>(x). Then, in the result of the standard replica calculations (which again inevitably involve the same 
cheeting as in the case of the previous model) we find that the corresponding free energy distribution function is not 
positively defined, which, of course, makes no physical sense. The point of this little methodological (and pedagogical) 
exercise is to demonstrate that it is not the replicas which are always responsible for all the troubles in the disordered 
systems world. In fact, the result of the replica calculations honestly reproduces the pathological nature of the 
original model itself: due to approximation, eq. (118p , we obtain the replica theory which does not correspond to any 
physical system. One can show that there exists no (positively defined) Gaussian distribution function of the random 
potentials V[<f>, x] which would provide the correlation function, eqs.((2]), in the "parabolic" form, eq. ([T8"l) ■ Moreover, 
the proper definition of the disorder potential is subject to important general constraints [22j regarding the shape of 
the correlation function U((j>), and neglection of these constraints may lead to unphysical results. 

Indeed, consider a random potential V(<p) and its Fourier representation V{p) = J dip V((f>) exp(— ip4>). Then the 
Gaussian distribution function of the random function V(p) has the form 

where the positive function G(p) is related with the correlation function U(<f>) via 

U{4>) = J |^G(p)exp(*p# (20) 
Expanding both sides of the above relation in powers of <f>, 

we can compare coefficients: the 2fc-th derivative of U (<j>) in the origin relates to the integral J dpG(p)p 2k which is a 
positive quantity. Hence, we have to be careful in our choice of the correlator U(<f>): if we truncate the expansion of 
U(<p) beyond some k*, such that C/( 2fc )(0) = for k > k* , we impose the condition 

J dpG(p)p 2k = for k > k*, (22) 

which cannot be satisfied for a positively defined G(p). Obviously, choosing the correlator U (</>), in the "parabolic" 
form, eq. (|18[) is in severe conflict with this constraint. If nevertheless, we would substitute such a parabolic correlator 
into the replica Hamiltonian and perform all the standard calculations, first we would see no apparent indications 
telling that something went wrong. Moreover, one can easily calculate the average free energy of this system to discover 
that it perfectly coincides with the one of the (physically consistent) random force model, eq. (fT7|) . The trouble appears 
in the calculation of more specific quantities. For example, the second cumulant of free energy fluctuations turns out 
to be negative which makes no physical sense. Including the next term oc <j> in the correlator's expansion can cure 
this problem, however, an inconsistency then shows up in the next order cumulants, etc. 
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The paper is organized as follows. In Section II we perform detailed analycis of the systematic replica field theory 
calculations for the random force model, cq. (fTT|) . and present the results for its free energy distribution function. In 
Section III similar solution is considered for the "parabolic" directed polymer problem described by the correlation 
function, eq. (|18p . providing not positively defined free energy distribution function. Finally, the subtleties of the 
replica method are discussed in Section IV. 

II. RANDOM FORCE MODEL 

Explicitly, the replica partition function, Eq.(|6]), of the system described by the Hamiltonian, Eq.(|17p. is 



+00 



Z(n,L) = JJ / dy a I T^^aix) exp 



4>a(0) = 



L n 



(23) 



Since the random parameters V\(x) and Vq(x) have Gaussian distribution with correlations defined in ea. (|15p . the 
disorder average (...) in the above equation is very simple: 



exp 



P dx^{Vi(x)4> a {x) + Vo(x)} = exp -0 1 \ dx V {u(f>a{x)<j) b {x) + v} 
J ° 0=1 J l Z Jo 0.6=1 



(24) 



Thus, the replica partition function, Eq. (|23[) . can be represented in the following form 



"> r-\-oo 

Z(n,L) = e^ 2 " 2 ^n / dy a y[y;L] 



where the "wave function' 1 



*[y;£] 



<t>a(L)=Va 



V(j) a (x) exp -j3H n [<f>] 



'0a(O)=O 

is defined by n-component scalar fields replica Hamiltonian 

(■L 



(25) 



(26) 



1 f I 

H n [4>] = - / dx\^2,[d x (j) a {x)\ - /3u^2 4>a{x)4>b{z 

Z Jt) \a = l a,b 

1 f L " 

= 2 / dX ^ ( f )a ^ Uab 

"' a,b=l 



ab 4>b(x) 

L J ° a,b=l 

The matrix 

U ab = d 2 x 5 ab + f3u 
can be easily diagonalized. It has (n — l)-degenerate eigenvalue 

Ai = d 2 x 

with (n — 1) orthonormal eigenvectors such that 

71 

J2ti=° (i=l,.-.,n-l) 

a=l 

and one non-degenerate eigenvalue 

A 2 = dl + /3nu 



(27) 



(28) 



(29) 



(30) 



(31) 
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with the eigenvector 



(32) 



Thus, the matrix U a b, Eg. (1281) . is diagonalized by the orthonormal transformation which is defined by the (n x n) 
matrix £? (such that X)"=i = % an< ^ £i=i = ^ab)- ln terms of new fields 



WW = 



(33) 



the Hamiltonian, Eq.(|27p. takes the form 

n-l 



1 /* 2 1 /* r 2 i 

#n [v] ""27 ^E^^W] + 2 i ^{[^^(a;)] -/3nM^(a;)| 



with the boundary conditions 



(34) 



(L) = nfy) = 



Q=l 



Correspondingly, the wave function, ea.(|26p. factorizes into 



*[y;£] 



II*o[ri(y);i] 



8=1 



*i[r„(y);i] 



where (with the proper choice of the integration measure) 



ip(L)=r 



Vo(r;L) = / T>tp(x) exp[--p dx[d x <p( 



<p(0)=0 



1 



2ttL 



/3 2 
exp( -— r 



and 



*i(r;i) 



<p(L)=r 



T><p(x) exp (-- / da;[/3(a x ^(a:)) 2 -/3 2 7iV^)] 



(35) 



'v(o)=o 

According to the above definition the wave function ^i(r, t) satisfies the imaginary-time Schrodinger equation 



with the initial condition 



*x(r;t = 0) = <5(r) 



(36) 



(37) 



(38) 



(39) 



(40) 



Eq. (l39j) describes the movement of a particle in the "reversed" parabolic potential. One can easily check that the 
solution of this equation satisfying the above initial condition is 



*,.(,-: I) : 1,(1) exp) ~a(t)r< 



where 



p_ (Ai 2 ) 1/4 
2irt 



(41) 



(42) 



and 



rr COs(V Xt 2 ) 

a(t) = /3VA — , ' 



(43) 
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where we have introduced the parameter 



A = j3nu 



(44) 



It should be stressed that the above solution exists provided 



o < xr < 



(45) 



This restriction indicates that at a given value of the parameter A the elastic string described by the partition function, 
eq.p8p (which contains negative mass!) goes to infinity at the finite time t c = it / (2^/~X). Coming back to the original 
random force problem one can reformulate the above restriction in the other way: for a given system size (" time" ) L 
the corresponding replica partition function Z n is defined only if 



< n c {L) 



Aj3uL 2 



(46) 



At bigger values of n the replica partition function is simply not defined (it is formally divergent). Taking into account 
that the quantity 7T 2 /(4/3uL 2 ) can be easily made less than one (e.g. by taking L sufficiently large), while the replica 
parameter n is still should be kept positive integer, we see that the above restriction becomes the fatal point for the 
whole replica scheme of the calculations. Of course, it is tempting to claim that since the replica n enter the partition 
function, eq. (!38p . in the form of the analytic parameter f3 2 nu, we can analytically continue it to arbitrary real values 
in the interval < n < n Cl eq. (|46p . But the problem is that the expression for the whole replica partition function, 
cqs. p5|) . ([55]) . requires that the replica parameter n must be still integer. In fact, the restriction, eq. (|i51) . reflects 
simple mathematical reality: the replica partition function Z n of the considered system is divergent at n > n c (L). 
In particular, at L > n/y/Afju, the replica partition function does no exist (divergent) for all integers n = 1, 2, .... In 
this situation (from the point of view of the common sense) the application of the replica method program for the 
system under consideration looks completely impossible. Nevertheless, in terms of the modus vivendi of the replica 
method this problem is overcome in a very simple way. It is in this point that the "doublethink" begins: wherever 
the replica n appears in a form of an analytic parameter it is immediately considered as real and belonging to the 
desired region (n — > in spin glasses, or < n < n c in the present case), while at the same time wherever n can 
not made non-integer (in the summations or in the products) it is still considered as an integer (note that similar, 
although slightly more sophisticated trick is used in the replica symmetry breaking construction). Thus, we continue 
our calculations just plainly assuming that the condition, eq. (!46l) . is satisfied. 
Substituting, the solutions eq. (|37p . (14T1) - (|43|) into eq. ([36|) , we get 



*[y;i] = 



2ttL 



n/2 



exp 



-1 

2L^ ' 



WT ) '71 y) 



/? VAL 2 cos(a 
2L 



(47) 



Using the relations, eqs. ([35| . (1301) and (I321) . and taking into account that the matrix £f is orthonormal, after some 
efforts in simple algebra we obtain the following result for the n-particle wave function, eq. (|26|) : 



*[y;£] = 



/ %/a! 2 / g 

sin(v / Al 2 ) \ 2ttL 



n/2 



exp 



n 



2L^ Va 



Er 2 

a=l 



2Ln 



/Al 2 costVAL 2 ) 




(48) 



where it assumed that (sicl) XL 2 = fJnuL 2 < 7r 2 /4 whatever the values of j3,u and L are. Substituting this result 
into eq. f|25[) and performing simple Gaussian integration (which, taking into account orthogonality of the matrix £f is 
easier to do in terms of the parameters n using expression, eq. (|47p ). for the replica partition function we finally get 
sufficiently simple result 



Z(n,L) 



I 



\ cos 



exp 



1 



/3 2 n 2 vL 



(49) 



Next, for further implementation of the general program of the reconstruction of the free energy distribution function 
as it was described in the Introduction, eas.([3|)- (|lip . let us introduce parameter 



w 



f3nuL 



(50) 
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which is confined in the interval < w < 7r 2 /4. In terms of this parameter the general relation between the replica 
partition function Z{n,L) and the free energy distribution function Pl{F) (cf. eq.©), 

/+oo 
dF P L {F) exp(-f3nF) (51) 
-oc 



takes the form 



where 



Z L (w) = / dFP L (F)exp( —F) (52) 



IV I 

Z L (w) = Z{— — -,L) = exp 

y /3uL 2 y/cOS{y/w) 



2u 2 L 3 



w 2 



(53) 



Rescaling the free energy, F — uL 2 f, for the probability distribution function Vl{I) of the random quantity / we 
get the following relation 

/+oo 
dfV L (f) exp(-wf) (54) 
-OO 

(where Vl(I) — uL 2 Pl(uL 2 /)). Next, performing the analytic continuation of the function Zl(w), eq. (|5"5)) , from 
the interval < w < 7r 2 /4 to the complex half-plain, Re{w} < 7r 2 /4, (which is unambiguous operation) the free energy 
distribution function 'Pl(Z) can be obtained via the inverse Laplace transform 



(55) 



/+20O J 
-too 2ttj 

where the integration goes over the contour parallel to the imaginary axes such that Re{w} < 7r 2 /4. 
In the thermodynamic limit, according to ea. (|53l) . 

lim Z L (w) = Z*(w) = 1 (56) 

Thus, according to eq. ([55|) . for the distribution function of the rescaled free energy fluctuations of the infinite system, 
V*(f) = linii-^oo Vh{f), w e obtain the following (universal) result [2lj : 

J-ioD lm vcos(vw) 

The overall form of this function is shown in Figure 1. It is interesting to note that this function is identically equal 
to zero at / > 0. This is easy to understand by making simple mathematical analysis of the integral in eg. ([57]). Indeed, 
since at / > the function [cos( v / Iy)]~ 1 / 2 exp(w/) quickly goes to zero at w — > — oo, the contour of integration in the 
complex plane can be safely shifted to — oo, which means that V*(f < 0) = 0. The fact that in the thermodynamic 
limit the upper bound for the free energy of the system described by the Hamiltonian, eq. (|17[) is equal to zero can 
also be explained in terms of simple physical arguments. First we note that the typical value of the random constant 
term J dxVo(x) scales as L 1 / 2 , which means that its contribution to the rescaled free energy / ~ F/L 2 scales as L~ 3 ^ 2 
and vanishes in the thermodynamic limit L — > oo. On the hand, since the contribution of the "trivial" configuration 
4>{x) = in the elastic and the random force terms of the Hamiltonian is equal to zero, any deviation from this 
configuration (due to the actions of the random force) can only reduce the energy. The asymptotic behavior of the 
function V*(f) in the limits / — > — oo and / — > —0 can be easily estimated by the saddle-point integration to yield: 

V4f -> -oo) ~ exp(- and V.(f -> -0) ~ exp^^) . 

At finite system size L, according to eqs.([53|) and ([55]) . the free energy distribution function is given by 



where 



r+lo ° dw 1 

VlU) = I , ■ \ — f exp e(L)w 2 + fw (58) 



2?T« y/cOS(y/w) 



-1 -0.8 -0.6 -0.4 -0.2 0.2 

FIG. 1: The thermodynamic limit free energy distribution function V*(f) 



The plot of this function for several values of the parameter e(L) is shown in Figure 2. At small values of e (large 
L) the function Vh{f) becomes close to the universal distribution function V*{f) shown in Figure 1, while at large e 
(when v u 2 L 3 ) 7 7 l(/) becomes almost Gaussian, as in this case the free energy of the system is dominated by the 
Gaussian random constant term J dxVo(x) of the Hamiltonian, eq . (fTTl) . 




-2 Jo i 2 



FIG. 2: Free energy distribution function P L (f), eq.HU, for: (a) e = 0.001, (b) e = 0.01, (c) e = 0.1, (d) e = 0.5, 

The result, eg. (1551) . constitutes the complete solution of the random force problem defined by the Hamiltonian, 
eq. (fT7f . Could we call this solution "exact"? One one hand, the explicit cheating in the the derivation of the replica 
partition function, cq. p5|) . makes the status of the obtained result rather indefinite. On the other hand, it should 
be stressed that this kind free handling with the integer/non-integer status of the rep lica parameter n is just the 
routine trick in all replica field theory calculations in disordered systems (see e.g. (MUM)- In other words, the results 
abtained in this way should be accompanied by the label " in the framework of the replica approach" . Fortunately, 
in this particular case, due to Gaussian nature of the considered model the result for the replica partition function, 
eq. p9")l . (where n £ (0, n c (L)) is a real parameter) can be confirmed by independent calculation without the use of 
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replicas (lg|, which allows to claim that the above result, eq. (|58|) . for the free energy distribution function is indeed 
exact. 



III. DIRECTED POLYMERS WITH "PARABOLIC" CORRELATIONS OF THE RANDOM 

POTENTIAL 



The replica partition function of directed polymers described by the Hamiltonian (fl]) with non-local parabolic 
correlations of the random potentials, Eq. (fTg)) , can also be represented in the form of eqs. (j2"5)) - (l2l)l) with the replica 
Hamiltonian 



H n [<f>] = - dx [^[dMx)] 2 + -puY\4>a{ 



x) - <p b {x) 



dx ^ 0a{x) Uab <f>b(x) 



a,b=l 



with the matrix 



Uab = (pi - PrmjSab + Pu 



(60) 



(61) 



Following the same route as in the case of the random force model, eqs. ([28l) - (|34[) . we note that the above matrix has 
(n — l)-degenerate eigenvalue Ai = d 2 — fJnu with (n — 1) orthonormal eigenvectors £f constrained by the condition 
J22=i & = (i = 1, ...,n — 1), and one non-degenerate eigenvalue A2 = d% with the eigenvector £° = l/^fn. In terms 
of the new fields (fii(x) = J22=i ^i^a, the replica Hamiltonian take the form (cf. eq.(j3U)) 



H n [tp] = 21 dx^2{[d x tpi(x)] - 0nu$(x)j + - / dx[d x 4> n (x)]' 

J „- 1 •/ 



(62) 



Similarly to the calculations of the previous section, eas (l551) - P5|) . for the replica partition function we get the following 
result (cf. eq.flUJl): 



Z(n,L) 



J cos(^\/ PnuL 2 ) 



(n-l) 



exp —(3 2 n 2 vL 



(63) 



where, as in the case of the random force model, it is assumed that finuL 2 < ir 2 /4. In terms of the parameter 
w = j3nuL 2 instead of eq. (|53l) we obtain 



Zz,p(w) = a/cos(v / w) a/cos(v / s) 



exp 



2u 2 L 3 



(64) 



and the probability distribution function of the (rescaled) free energy fluctuation Vl(I) is given by the inverse Laplace 
transform (cf. cq . (f55j) ) : 



r+loc dw ~ 
Pl .(/)= / gUzH ex P (u,/), 



(65) 



Performing (numerical) integration in the above equation one can easily find that this function is not positively defined 
for any values of the parameters e = v / (2u 2 L 3 ) and k = 1/ (2/3uL 2 ). For example, in the zero temperature limit (when 
k — y 0) the above equation reduces to (cf. eq . (1551) ) : 



lim vlm = nu) = 



dw 
2~Tri 



cos(y/w) exp e(L)w 2 + fi 



(66) 



The plot of this function for several values of the parameter e(L) is shown in Figure 3. We see that unlike Vl(/) 
(Figure 2) of the random force model, the function V^(f) even at small L (large e) (when it is almost Gaussian) has 
always a kind of the negalive "kink" at sufficiently large /. 
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Unfortunately we are not able to confirm (or reject) the result, eq. (|63|) . by independent calculations, as the replica 
theory defined by the Hamiltonian, eq. {J5DJ), does not corresponds to any physical system (see the discussion of this 
issue in the Introduction, eqs. ([T9|) -(f22 j) ). On the other hand, the clear lesson which we can learn from the exercise 
considered in this Section, is that unlike the "honest cheating" with the status of the replica parameter n (discussed 
in previous Section), any approximations made at the stage of the replica calculations (the exact correlator U((f>), 
eq. (J2J) , is replaced by its truncated expansion, eq.(fT8j)) could be just fatal for the physical meaning of the obtained 
results. 




FIG. 3: Free energy "distribution function" P£(/) as it given by eq.i[66f for: (a) e = 1.5, (b) e = 0.5, (c) e = 0.05, (d) 
e = 0.0005, 



IV. DISCUSSION 

The standard program of the replica method is formulated as follows: first, for an arbitrary positive integer n we 
have to calculate the disorder average of the n-th power of the partition function, Z n = Z(n) which is expected to be 
an analytic function of the replica parameter n; second, we have to perform an analytic continuation of this function 
from integer for arbitrary real or complex values of n; and third, we have to take the limit n — > (if we are interested 
in the average free energy only) or we have to perform an integration over complex n (if we are deriving the free energy 
distribution function). This third step is usually accompanied by taking the thermodynamic limit, which assumes 
that the system size L is taken to infinity. The prescription of the replica method indicates that the two limits, n — s- 
and L — > oo, has to be taken simultaneously such that the product nL u (where an exponent w defines the scaling of 
the free energy with the system size) is kept finite. 

In fact, the whole experience of the replica calculations in disordered systems shows that except for trivial cases this 
program, as it is formulated above, is never followed! The typical illustration of the replicas realpolitik is provided by 
the studies in the mean- field spin glasses 0] ■ First of all, since the system is sufficiently complicated the computation 
of the replica partition function Z(n) can be done here only in the saddle-point approximation and not exactly. This 
makes further analytic continuation to non-integer n somewhat doubtful because the neglected terms which are small 
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at integer n and large system size could become essential in the limit n — > 0. Moreover, it turns out that at large 
n the replica partition function growth as ~ exp(n 2 ) which means that its analytic continuation to non- integer n is 
ambiguous. If nevertheless we would just plainly take the limit n — > in the obtained expression for Z(n) , we would 
get the so called replica symmetric (RS) solution which at low temperatures is unphysical since it reveals negative 
entropy and many other bad things. In view of the remarks made above, of course, this is not surprising. The strategy, 
which is called the replica symmetry breaking (RSB) scheme Q, and which is generally believed to provide correct 
results, is essentially different. In this scheme, all the above three steps, (computing Z(n), analytic continuation in n 
and the limits n — > and L — > oo) are performed simultaneously! Similar (although slightly simplified) scheme works 
perfectly well also in the case of the Random Energy Model of spin glasses @. For both models the replica results 
are confirmed by independent mathematically rigorous calculations [9(. There are many other systems for which the 
results of the replica calcualtions, although can not be confirmed rigorously, are generally accepted to be correct [23| . 

To understand what is going on sometimes it is useful to consider an example of very simple system. Replica 
calculations for the model studied in this paper does not involve any kind of the RSB " magic" . The model is so simple 
that initially one gets an illusion that every step of the calculations could be under good control. For instance, unlike 
the above example, the replica partition function here can be computed exactly for any (finite) system size, eq. (|49|) . 
Nevertheless, proceeding with further steps of the replica method program one finds that either the attempt should 
be aborted, or one has to start cheating again. Indeed, the result, eq. (|49[) . for the replica partition function turns out 
to be valid only for finite number of the integer points: n < n c = [ir 2 / (AftuL 2 )], cq. (|46[) , since at n > n c the quantity 
Z(n) is not defined (it is formally divergent). Moreover, for sufficiently large system size, L > L c (n) — n / (2y/ f3u) , 
the replica partition function Z(n) is not defined for all positive integer n (including n = 1). In other words, in this 
situation the replica partition function of the considered system, defined as Z n (where n is a positive integer) simply 
does not exist! But even if L < L c (n), so that Z(n) is still defined at finite number of the integer points, it is evident 
that its analytic continuation for non-integer values of n is completely ambiguous. If nevertheless, one neglects all the 
above observations and accept the result, eg. (1491) . as valid for all real n in the interval < n < n c , then everything 
becomes just fine. The analytic continuation of the function, eg. (14U1) . from the finite interval < n < n c to the the 
complex half-plain, Re{n} < n c is unambiguous, and in this way one obtains beautiful and correct results for the free 
energy distribution function (Figures 1 and 2). In this particular case one can be sure that obtained results are indeed 
correct as for the system under consideration the quantity Z n can be computed directly for any real n £ [0, n c ] [la ]. 

All the experience of last decades convincingly demonstrate that with a few exceptions the replica method does 
give correct results, and this can not be explained by simple coincidences. We know very well how it works, and we 
do know that the replica calculations inevitably involves cheating. The question is then, why it works? 
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